Two-nucleon correlation effects in knockout reactions from ^^C 
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Reactions that involve the direct and sudden removal of a pair of like or unlike nucleons from a fast 
projectile beam by a light target nucleus are considered. Specifically, we study the three two-nucleon 
removal channels from ^^C that populate final states in the ^"Be, and ^"C reaction residues. The 
calculated two-nucleon removal cross sections and the residue momentum distributions are compared 
with available high energy data at 250, 1050 and 2010 MeV per nucleon, data that are inclusive with 
respect to the bound final-states of the residues. The measured np-removal cross sections only are 
significantly greater than the values calculated, suggesting that the reaction mechanism observes 
enhanced np spatial correlations compared to those present in the shell-model wave functions. 
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I. INTRODUCTION 

There is an extensive literature on the development 
and the use of high-energy electron-induced knockout of 
both one or two nucleons from the ground states of stable 
nuclei as probes of nucleonic correlations and of pair cor- 
relations. Relevant and recent reviews include references 
[ll-U , which cite the essential experimental developments 
and data sets: see also the citations in [F]. Quite recently, 
new conclusions have been drawn regarding the obser- 
vation of enhanced short-ranged correlations (SRC) be- 
tween nucleon-nucleon pairs from (e, e'pN) experiments 
on a carbon target ^J. Specifically, measurements that 
selected high momentum transfer and large missing mo- 
mentum events from the few-body final-state phase-space 
found evidence that np pairs were more than an order 
of magnitude more prevalent than like-nucleon pairs in 
this selection; deduced from the corresponding (e, e'np) 
and (e, e'2p) yields. The probabilities of these high rela- 
tive momentum two-nucleon components have been com- 
puted, using the microscopic ground-state wave functions 
from ab-initio variational Monte Carlo calculations, for 
systems with masses A < 8 [aj7| and, using a linked 
cluster expansion, for A > 12 Q. These computations 
involve the wave functions of the nucleons over the entire 
volume of the nuclei concerned. These theoretical studies 
are strongly suggestive that the tensor force, in the spin 
S — 1, isospin T = 0, np channel, plays the major role 
in generating these enhanced np-spatial correlations and 
their associated high-relative-momenta signatures in the 
two-nucleon density distribution. In such (e, e'pN) mea- 
surements, the specific final-state phase space selection 
(of high relative momentum back-to-back two-nucleon 
events) provides the leverage and the amplification of the 
short-ranged pair-correlation sensitivity. 

The behavior of the wave functions of single nucleons 
and of pairs of nucleons (in a mass A + 2 projectile) can 
also be probed if the nucleons are removed (suddenly) 
in fast collisions with a light target nucleus. Within this 
strong-interaction probe case, the sensitivity is now to 



the nuclear wave functions at and near the nuclear sur- 
face. Such processes are also referred to as one- and 
two-nucleon knockout (or removal) reactions. These di- 
rect reaction mechanisms, combined with 7-decay spec- 
troscopy, to determine the final state of the reaction 
residues, and their momentum distributions, are being 
actively exploited as a spectroscopic tool to study the 
evolution of nucleonic single-particle structure near the 
two (often very displaced) Fermi surfaces of exotic nuclei. 
They are proving to be robust techniques, see e.g. 

Unlike the final-state-exclusive measurements in elec- 
tron induced reactions, observables in such nuclear- 
induced fast nucleon removal reactions are inclusive with 
respect to the final states of the removed nucleon(s) and 
of the fate of the struck light target nucleus. The cross 
sections are thus relatively large. Measurements usu- 
ally consist of the total removal reaction yield, often 
the momenta of the fast, forward travelling projectile- 
like residues, and sometimes the differential yields to 
the ground and bound excited states of the mass A + 1 
or mass A residual nuclei; the latter obtained by 7-ray 
spectroscopy. However, in such collisions between the 
projectile and a light composite target nucleus, such as 
beryllium or carbon, the reaction is geometrically very 
selective [13], and the removal of two nucleons will be 
enhanced if nucleon pairs have a strong spatial corre- 
lation (and localization) in the projectile ground state. 
An interesting question therefore is whether like and un- 
like two-nucleon removal under such conditions also ex- 
hibits any evidence of enhanced np over nn and pp spa- 
tial correlations on a longer length scale than implied by 
the (SRC) observations of the electron knockout data? 
Here, the spatial (geometrical) selectivity of the reaction 
mechanism would provide the leverage and probe of the 
presence of spatially localized pairs and evidence of an 
enhanced np correlation. 

Experimental data are available in the form of high- 
energy, primary beam measurements of the inclusive 
cross sections to the bound states of the residues after np, 
nn and pp removal [l^, [Si • We will show that these data 
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reveal a significant enhancement of unlike-pair yields, 
(J-np, relative to those, (T_n„ and (T_pp, for like-nucleon 
pairs, the enhancement being significantly greater than 
would be expected (trivially) from the numbers of such 
pair combinations available. 

The removal of two (well-bound) nucleons of the defi- 
cient species from asymmetric nuclei has been shown to 
proceed as a direct reaction [l^l- A description of the 
reaction making use of configuration-mixed shell-model 
wave functions and elastic and inelastic breakup contri- 
butions to the removal cross section has been shown to re- 
produce experiment Such calculations have also re- 
cently been extended to describe the reaction-residue mo- 
mentum distributions [l6l [T7| , where it was shown that 
the shapes and widths of the residue momentum distribu- 
tions are indicative of both the total angular momentum 
(/) and the total orbital angular momentum components 
{L = ti + I2) oi the removed nucleon pair 

Here we will exploit the eikonal reaction model in the 
isospin formalism [l^, [13 for the removal of like 
(T = 1) and unlike (T = 0, 1) nucleon pairs. We also dis- 
cuss the possible role of indirect population of the residue 
final states of interest, i.e. by single-nucleon knockout 
paths that populate particle-unbound states of the inter- 
mediate mass A-\-l system. Such indirect paths severely 
limit the number of systems that, in the absence of an 
empirical means to distinguish between direct and indi- 
rect knockout events, may be studied quantitatively with 
unlike pair knockout. These indirect contributions to fi- 
nal state yields are expected, in general, to be large in 
np removal from asymmetric systems [20| where one or 
other of the nucleon thresholds will be at relatively low 
excitation energy. 

Attractive np-removal test cases thus suggest symmet- 
ric light nuclei, e.g. ^^C or ^^O. For such less-massive 
cases one must however consider core-recoil effects in the 
reaction and the analogous center-of-mass corrections to 
the (fixed center) shell-model two-nucleon amplitudes. 
Both are expected to affect, in detail, the absolute mag- 
nitudes of the cross sections. These are discussed here. 
Regardless of these two (relatively small) effects, compar- 
isons of the relative like- and unlike-pair knockout yields 
may be made. 

In Section we outline the specific features of the 
reactions using ^^C projectiles. The necessary formalism 
has been presented elsewhere and will only be outlined 
in Section IIIIl exploiting the notation used in previous 
work. The calculated results for ^^C are discussed in 
Section [IVI and a summary made in Section Ivl 



II. CARBON INDUCED REACTIONS 

Our current expectation is that a quantitative discus- 
sion of pair-correlation effects can only be made when 
the residue final states are populated predominantly via 
a single-step direct reaction. It is thus vital to select ex- 
amples that minimize the indirect (evaporative) contri- 



butions to the reaction yield: and that we assume cannot 
(at present) be distinguished from the direct reactions of 
interest. We must therefore consider systems with large 
(and symmetric) nucleon separation thresholds, such that 
the single-nucleon removal strength to particle-unbound 
states will be weak. For the same reason, we require the 
projectile nucleus to be relatively light to minimize popu- 
lation of the (bound) final residue states via the removal 
of deeply bound (non-valence) nucleons. Ideally, any ex- 
periment would also measure single-nucleon and like-pair 
removal in addition to the unlike pair removal, in order 
to verify shell- and reaction-model predictions for the dis- 
tribution of single-nucleon spectroscopic strength and for 
the direct T = 1 pair knockout. The above considerations 
severely constrain the potential candidates, the best ex- 
amples being ^^C and ^^O. The former of these will be 
considered in detail here. The choice of an = Z pro- 
jectile poses the additional (experimental) complication 
of distinguishing the mass A residue from the incident 
beam, the residue and projectile having identical mass- 
to-charge ratios. 

Consideration of two-nucleon knockout from ^^C is 
valuable for two reasons. First, its shell model descrip- 
tion and that of the residual nuclei ^°C, ^°Be and 
are extensively studied and so establish a valuable point 
of reference. Second, the existing experimental cross sec- 
tions for two-nucleon knockout from ^^C [H, are ac- 
curate to w 10% and were taken at high energies, where 
the eikonal model used here is at its most reliable. The 
experimental data were obtained using reactions of a car- 
bon beam on a carbon target at 250, 1050 and 2100 McV 
per nucleon incident energies. These show an expected 
enhancement of the ^"^B production cross sections (np re- 
moval) over those for the nn (to ^"C) and pp (to ^°Be) 
removal reactions (see Table Hill . The significance of this 
observed enhancement is quantified here. 

In addition to two-nucleon removal, single-nucleon re- 
moval cross sections were also measured at high energy 
(and previously studied theoretically, see Ref. 10]). 
These verify the expectation, from the shell-model, that 
the majority of single-particle removal strength is ex- 
hausted in transitions to final states lying below the nu- 
cleon separation thresholds of ^^C and ^"'^B. This point 
will be developed further in later Sections. The relevant 
one- and two-nucleon (and a-particle) separation thresh- 
olds are illustrated in Fig. [TJ 

The simplest of the three two-nucleon removal cases 
is two-neutron (nn) removal, leading to the resid- 
ual nucleus. This has only two bound final states below 
the first proton threshold, the 0+ ground state and a 2+ 
excited state at 3.354 MeV [21]. The level scheme of 
^*^Be, the residual nucleus in the pp removal case, is only 
slightly more complex. We must consider population of 
the 0+ ground state, the 2+ states at 3.368 and 5.958 
MeV, and a second 0+ state at 6.179 MeV, aU below 
the neutron threshold of 6.812 MeV. The most compli- 
cated final state is that for ^°B, the np knockout residue. 
The p-shell shell-model calculations used include even- 
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FIG. 1: One- and two-nucleon and a-particle separation 
thresholds relevant to nucleon-removal reactions from ^^C. 
The '^^C two-nucleon separation energies are S2p ~ 27.184, 
S2n = 31.184 and S„p = 27.412 MeV. Proton evaporation 
from the single-neutron removal residue ^^C would be ex- 
pected to be the largest indirect pathway, but the shell-model 
calculations suggest very little spectroscopic strength to states 
above 8.689 MeV in "C. 



parity states up to a maximum of spin 3. For ^"B we 
have shown states up to the proton separation threshold, 
however the low a-particle separation threshold means 
that the T = states at 4.774 (3+), 5.180 (1+) and 5.920 
MeV (2+) are reported to decay (with branching ratios 
of near 100%) by a-emission. We will show the cross sec- 
tions for population of these states but will assume that 
they do not contribute to the calculated yield. The 
T = 1, 2+(5.184 MeV) state is reported to have a 16% a 
emission branch. These "'^^B states are illustrated in Fig. 
[5J We note that the ^"B spectrum also contains several 
negative parity states which are not expected to be pop- 
ulated by the nucleon removal reaction mechanism. 

Direct, unlike-pair-removal cross sections are expected 
to be larger than those for like-pairs. Assuming the sim- 
plest p— shell 7r[0p3/2]^i^[0p3/2]^ structure for ^^C, the pn- 
format TNAs for pair removal to a residue final state 
with spin Jf are given by the appropriate coefficients of 
fractional parentage, with value y^2Jy~+T. The cross 
sections to a given final state, being proportional to the 
square of these TNA, suggest unlikedike two-nucleon in- 
clusive cross section ratios cr-np/<J^2N = 16/6 « 2.7 if aU 
strength leads to bound final states. The same result is 
obtained from the counting of available valence-nucleon 
pairs. 

The primary motivation for the present study is that 
this simple combinatorics expectation is not corroborated 
by the available data; the ratio anp/cpp {(Jnp/<Jnn) from 
the experiments at 250, 1050 and 2100 MeV/nucleon are 
8.1 (8.9), 5.3 (6.3) and 6.0 (8.5), respectively [H, the 
observed enhancement of the np removal cross section 
being significantly larger than expected. 
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FIG. 2: States of the mass A — 10 residues populated in the 
two-nucleon knockout. The spin and isospin labels {Jf,Tf) 
are indicated. All states included are of positive parity. Lev- 
els assumed to be part of the isospin multiplet are connected 
by dashed lines. The lowest particle thresholds are also in- 
dicated. States above the a-particle threshold in ^''B are ex- 
pected to decay via a emission, with the exception of the 
5.164 MeV, T = 1, J'" = 2+ state which has an 84% 7-decay 
branch. 

Clearly our simple estimate, from the number of avail- 
able pairs of each type, assumes that all of the removal 
strength leads to bound configurations in the residue of 
interest, which is not the case, and part of the enhance- 
ment noted above may be attributable to this (channel- 
dependent) fraction of events leading to particle-unbound 
states. Here we seek to address this question quan- 
titatively, by use of a direct reaction model for two- 
nucleon removal. Specifically, our aims are to quantify: 
(i) whether this enhancement is accounted for by those 
two-nucleon correlations that are included in the trun- 
cated, p-space shell model calculation, and (ii) to make 
a simple estimate of the contributions one might reason- 
ably expect from indirect reaction pathways in this model 
picture. 

III. FORMALISM 
A. Two-nucleon overlap 

The formalism used is based on that developed in Refs. 
[TtI . Isospin- format two-nucleon amplitudes (TNAs) 
will be used and are expected to offer a good description 
of the light symmetric systems considered here. Thus, in 
the unlike pair removal case, we assume a common set of 
nucleon orbital wave functions, to be discussed below. 

We evaluate the cross sections for transitions from the 
projectile initial (ground) state i, with spin {Ji,Mi), to 
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particular residue final states /. The residue is assumed 
to be a spectator in the sudden reaction description and 
its state is not coupled to the reaction dynamics. The 
direct reaction will then probe the two-nucleon overlap. 
We denote the A-body final states by where the 

label F = (f,Mf) includes the angular momentum pro- 
jection Mf. The two-nucleon wave function (two-nucleon 
overlap) of the removed nucleons 1 and 2 is written as 
^'^^^(1,2), where 

—^JiMiT,T,\' ' 

^($(^)(A)|vl/,(A,l,2)) 
= ^ &J{I^JiJfMf\.hM,) 

Ifj,Ta 

{TrTfTf^n) [J^M®^hM]l; ■ (1) 
I 




FIG. 3: (Color online) Diagram of the angular momentum and 
isospin couplings, and the coupling orders assumed. Specifi- 
cally, the LS-couplings used in this paper are presented. 



with Da = l/-\/2(l + (J/Si/Ja)- The angular momentum 
and isospin couplings used are summarized in Fig. [31 
The nucleon-wave functions (l>^j{i) are 

c^^^.^ii) = U0[n)Ytx{h) ■ (3) 

We note that no explicit account is taken of additional 
(short, medium, long range) correlations other than those 
contained within the shell-model description, and that 
enter Eq. ([T]). These are seen to arise from: (a) antisym- 
metry and angular momentum coupling of the nucleon 
pair, and (b) the (shell-model) two-nucleon overlap, via 
the weights and phases of the contributing TNAs. If 



Here the set of available active two-nucleon configura- 
tions, with counter a, consists of particular pairs of 
single- particle orbitals [/3i, /32], where the index /3 = 
{nHj) denotes the state's spherical quantum numbers. 
The C^"^ are the two-nucleon amplitudes (TNAs) that 
express the parentage of the residue final state, when 
coupled to a particular two-nucleon configuration a, in 
the projectile ground state. Here, as elsewhere, these are 
taken from a truncated-basis shell model calculation. 



Expressed in iS'-coupling, the antisymmetrized two- 
nucleon wave function of Eq. ((1} is 



(2) 

I 

there are significant additional strong-interaction gener- 
ated pair correlations, missing from our description, then 
we might expect empirical cross sections to deviate sig- 
nificantly (and be enhanced) relative to the shell-model- 
correlated model used here. 



B. Centre-of-mass corrections to the TNA 

Spectroscopic factors calculated within the fixed-center 
shell-model basis require a center-of-mass correction fac- 
tor to be applied [24I. Essentially, the shell-model two- 
nucleon overlaps are calculated relative to the centre of 
mass of all A -I- 2 nucleons and not relative to the re- 
maining A nucleons. Such corrections were discussed by 
Pinkston (23l - l25| (see also Ref. [1^). In previous, pre- 
dominantly sd-shell applications such corrections were 
expected to be small, but for the examples discussed here 
they may be more significant, particularly when one has 
an ambition to compare absolute cross sections with ex- 
periment. 

In the simplest (harmonic oscillator) limit, the require- 
ment is to multiply the shell-model TNAs by a factor 
[{A + 2)/A]^/^, with N the number of oscillator quanta 
of the orbitals of the active nucleons. In the present 
case, of "'^^C(— 2N), this would enhance all cross sec- 
tions by ~ 20%, independent of the reaction channel 
and the residue. Thus, conclusions regarding the relative 
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strengths of the different two-nucleon removal channels 
will be unaffected and we do not apply these corrections. 
Clearly, a more complete understanding will be necessary 
for a precise consideration of the absolute cross sections. 

C. Reaction dynamics and approximations 

We exploit eikonal reaction dynamics with the required 
nucleon- and residue-target S'-matrices calculated in the 
optical limit of Glauber theory |27i, i28j. These elastic 
S'-matrices are thus calculated assuming that the residue 
and nucleons travel on straight line paths through the 
interaction field of the target. The residue-target S'- 
matrices are calculated by double-folding their densities 
with an effective nucleon-nucleon (NN) interaction. The 
absorptive nature of these interactions naturally local- 
izes the reaction to the projectile surface flT\ . A detailed 
derivation of the two-nucleon removal cross section is pre- 
sented elsewhere as is the formalism for their mo- 
mentum distributions [13, [13 ■ These are not reproduced 
here. 

The key approximations are as follows. The removal 
of nucleons is sudden, i.e. their co-ordinates are assumed 
to vary slowly compared to the timescale of the reaction, 
and can be assumed frozen during the interaction. In 
previous work we have also (reasonably) made the no- 
recoil (heavy residue) approximation in the diffraction- 
stripping terms of the removal cross section [Tsj ; i.e. we 
have assumed coincidence of the residue and projectile 
centers of mass and impact parameters, be ~ b. These 
core recoil effects will certainly introduce corrections for 
the lighter residues in this work, however, as will be quan- 
tified later, at the high projectile energies considered here 
these diffractive-stripping terms contribute a significantly 
smaller fraction of the cross section than the two-nucleon 
stripping terms. Thus, recoil will have a relatively minor 
effect on the calculated removal cross sections and the 
discussion of this paper. This effect will be fully quanti- 
fied in future work. 



IV. NUCLEON KNOCKOUT FROM ^^C 
A. Methodology 

The experimental information available consists of fi- 
nal state inclusive cross sections to bound states of the 
mass A = II, 10 residues and some associated residue 
momentum distribution information. We are able to 
identify each known positive parity final state (A = 10) 
and negative parity final state {A = II) with a well- 
defined shell model state. The shell-model calculations 
were performed with OXBASH [2^. We will consider the 
results when using the WBP interaction 1301 . and also 
compare these with the pjt interaction [31[ outcomes, 
in the p— shell model space. The former interaction was 
used in previous work on nucleon removal from ^^C [lo| . 



The shell- model spectroscopic factors C'^S and the two- 
nucleon spectroscopic amplitudes, the C^^ , as enter Eq. 
dD), were calculated for each residue final state. To dis- 
cuss possible indirect processes we must consider the pop- 
ulation of highly-excited, unbound states in the single- 
nucleon removal cases. 

The p— shell model space allows removal of nucleons 
from the active 2 andpi/2 single-particle orbitals. The 
radial wave functions used to describe these nucleonic 
states are calculated in a Woods-Saxon potential well. 
The radius and diffuseness parameters were fixed to be 
To = 1.310 fm and a = 0.55 fm, as were used in the single- 
nucleon knockout analysis of Ref. In addition, we 

include a spin-orbit potential of depth Vso = 6 MeV. The 
ground-state to ground-state nn, pp and np separation 
energies, S2N, from ^^C are 31.841 MeV, 27.184 MeV and 
27.41 MeV, respectively. The nucleon bound state wave 
functions for each reaction and transition, with final-state 
excitation energy Ef, were calculated using effective nu- 
cleon separation energies of 

Bf^SN + Ef, (4) 

for single nucleon removal and 

Bf = {S2N + Ef)/2 (5) 

for each nucleon in two-nucleon removal. In the np- 
removal case we assumed a common set of radial wave 
functions for the neutrons and protons. These were cal- 
culated as above but assuming each nucleon carried a 
charge of 0.5e. 

The required nucleon- and residue-target S'-matrices 
were calculated using the optical limit of Glauber theory 
(27} - i.e. by double folding each constituent particle and 
target density with an effective nucleon-nucleon (NN) 
interaction. The essential input parameters, given the 
dominance of the stripping mechanism, are the free nn 
(= pp) and np cross sections, taken from the parametriza- 
tion of Ref. [13], and the residue and target densities. 
The former are given in Table HI The effective NN in- 
teraction was assumed to be zero-range, represented by 
a (5-function, and the reakimaginary forward scattering 
amplitude ratios a„„ and a„p were obtained from a poly- 
nomial fit to the values tabulated in Ref. [1^. These 
are listed in Table U for the three energies of interest. 
The point-nucleon density distributions of the target and 
the residues were assumed to have Gaussian shapes, with 
root mean squared (rms) radii consistent with Glauber- 
model analyses of the measured interaction cross sections 
[S^l- For simplicity, we assume a single matter rms ra- 
dius for each mass nucleus; for ^^C target and A = 11 
and A = 10 residues, the mass radii were 2.32, 2.11 and 
2.30 fm respectively. The latter value is that quoted for 
i"Be in Ref. [H. 
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TABLE I: Nucleon-nucleon efTective interaction parameters 
used for the calculation of the nucleon- and residue-target 
S-matrices. 



Energy (MeV/u) 


cr„„ (mb) 


cr„p (mb) 






250 


22.3 


36.6 


0.94 


0.49 


1050 


49.6 


41.6 


-0.08 


-0.46 


2100 


63.9 


46.5 


-0.20 


-0.47 



B. Single-nucleon removal cross sections 

We first briefly review the single-nucleon removal re- 
sults, as were previously discussed in Ref. Of par- 
ticular interest here is the proportion of spectroscopic 
strength exhausted below particle separation thresholds, 
and hence the shell-model prediction of single-particle 
strength lying above the nucleon separation thresholds. 
For the wbp interaction 3.93 units of spectroscopic 
strength are associated with states below these thresh- 
olds with the remaining 0.07 fragmented over many 
states above 10 MeV in excitation. Calculations using 
the PJT interaction distribute the single-particle spectro- 
scopic strength similarly, with 3.97 units associated with 
states below the mass A = 11 nucleon thresholds. 

The results for single-nucleon removal are very similar 
to those presented previously [13] ■ The small differences 
can be traced to our inclusion of a spin-orbit term in 
the nucleon bound states potential. The detailed decom- 
position of the single-particle cross sections for each fi- 
nal state into their stripping and diffraction components, 
for the beam energy of 2100 MeV/nucleon, are shown in 
Table [TTl A discussion of the single-nucleon shell-model 
suppression factors can be found in Ref. [l3|. Here we 
simply note that the theoretical cross section overesti- 
mates the experiment by approximately a factor of two, 
consistent with the results from electron-induced proton 
knockout [1]. 

With regards to possible indirect contributions to np 
knockout, we can estimate an upper bound assuming that 
the remaining shell-model C^S strength (of 0.07) from 
both proton and neutron channels will populate the np 
residue indirectly. Doing so we obtain an indirect contri- 
bution of w 3 mb. It is however unlikely that this entire 
residual p-shell strength would lead to the ^"B residue. 
The removal of deeply-bound Osi/2 nucleons may also 
lead to the ^"B residue, though the single-particle cross 
section will be significantly smaller than those presented 
in Table nil due to the stronger binding. We return to this 
indirect contribution discussion in the following Section. 



C. Two- nucleon removal cross sections 

The stripping mechanism component of the two-like- 
nucleon removal cross sections were previously discussed 
in Ref. (Sal . Here we extend these results to consider the 



TABLE 11: Single-nucleon knockout cross sections (mb) for 
the ^^C projectile incident on a carbon target at 2100 MeV 
per nucleon. The single-nucleon removal cross sections pre- 
sented, (J-iAT, include both the spectroscopic C^S (from the 
WBP interaction) and center of mass correction factors (here 
12/11). The results are very similar to those of Ref. [l^. The 
inclusion here of the nucleon bound states spin-orbit potential 
accounts for the very minor differences observed. 



Residue 




<^ str 




C^S 


U-IN 


"C 


3/2" 


20.50 


2.17 


3.16 


78.18 




1/2- 


18.45 


1.81 


0.58 


12.82 




3/2" 


18.74 


1.81 


0.19 


4.26 










sum 


95.36 










exp. 


46.50 ±2.30 


"B 


3/2" 


21.11 


2.30 


3.16 


80.70 




1/2- 


18.86 


1.90 


0.58 


13.14 




3/2" 


19.09 


1.88 


0.19 


4.35 










sum 


98.81 










exp. 


53.80 ±2.70 



stripping-diffraction mechanism contributions and con- 
sideration of unlike-pair knockout to confront the avail- 
able experimental data quantitatively. As was discussed 
earlier, disregarding the configuration mixing inherent in 
the shell-model two-nucleon overlap, we might expect the 
unlike pair knockout yields to be enhanced by approxi- 
mately a factor of 16/6 over those for like nucleon knock- 
out, from the counting of available [^3/2]^ pairs. As was 
also pointed out, this zeroth-order estimate assumes all of 
the two-nucleon strength will fall below the first particle- 
emission thresholds. In practice the number of bound 
states and the fraction of the removal strength leading 
to bound states is different for the three residues. The 
lowest relevant particle separation thresholds are as fol- 
lows: i°C, Sp = 4.006 MeV; i°Be 5'„ = 6.586 MeV; 
^°B, Sa = 4.461 MeV. Both the number of pair com- 
binations and the distribution of two-particle removal 
strength should be reasonably accounted for when using 
the shell model two-nucleon amplitudes. 

Results for the inclusive cross sections for two-nucleon 
removal from ^^C to the three residues, at beam ener- 
gies of 250, 1050 and 2100 MeV per nucleon, are given 
in Table IIIII The final-state-exclusive like and unlike 
two-nucleon removal cross sections, and their decompo- 
sition with the contributing reaction mechanisms (i.e. 
stripping, a str, diffraction-stripping, Uds, and estimated 
diffraction, (idif) are shown in Table ITVl for the 2100 MeV 
per nucleon case and when using the wbp interaction 
TNAs. Table IIIII shows the calculated theoretical like- 
pair removal cross sections atin to and ^''Be, are in 
reasonable agreement with the experimental data, Uexpi 
of Refs. 

From Table ITVl we note that, at this projectile energy, 
the two-nucleon stripping (absorption) term, Ustr, 
counts for «70% of the calculated cross section. The 
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TABLE III: Theoretical and experimental cross sections for two nucleon knockout from ^ C, for projectile energies of 250, 1050 
and 2100 MeV per nucleon. All cross sections are in mb. The TNAs used were calculated using the wbp interaction. 



Energy 




i»Be 






10(j 










MeV/u 


O'th 




O'exp/o'th 


O'th 


^exp 


Oexp 1 Oth 


Oth 


(^exp 


Oexp/ Oth 


250 [12] 


7.48 


5.88±9.70 


0.79±1.30 


5.80 


5.33±0.81 


0.92±0.14 


21.57 


47.50±2.42 


2.20±0.11 


1050 [13| 


6.62 


5.30±0.30 


0.80±0.05 


5.13 


4.44±0.24 


0.87±0.05 


19.27 


27.90±2.20 


1.45±0.11 


2100 [13] 


6.52 


5.81±0.29 


0.89±0.04 


5.04 


4.11±0.22 


0.82±0.04 


19.02 


35.10±3.40 


1.84±0.18 



cross sections for those T = 1 states common to all three 
residues, namely the first 0+ and 2+ states, are also es- 
sentially equal, the minor differences in the calculations 
arising from the small differences in the separation ener- 
gies for each system. 

For np removal to ^''B, the cross sections are shown 
for the nine p-shell shell-model final states below the first 
nucleon threshold. However, the first 2+, T = 1 state is 
known to decay by a-emission with a branch of = 16%. 
For this state, this branching has been accounted for in 
the a-2N value presented. The cross sections are also 
shown for the three highest energy T = states, but 
these states are reported to decay (with — 100%) by 
a-emission. These cross sections are not included in the 
summed ^"B yield. These exclusive ^°B cross sections 
will be discussed further in Section flV DI and in Figure [Sj 
in connection with their associated momentum distribu- 
tions. We note here that if one uses the pjt shell-model 
interaction TNAs then the summed ^°B yield is 18.73 
mb, compared to the 19.02 mb shown in Table IIVI The 
calculated inclusive cross sections are thus robustly de- 
termined within the p-shell model space calculations. 

As for the present like-nucleon cases, previous analyses 
of like two-nucleon removal data for asymmetric sd— shell 
nuclei have indicated a suppression of the measured val- 
ues compared to the theoretical model being used. Previ- 
ously this factor was of order Rs{2N) — (Texp/crth ~ 0.5. 
The present higher energy, like-nucleon calculations sug- 
gest a larger value of Rs , see Table IIIII As was noted 
earlier, the present calculations do not include center-of- 
mass (recoil) corrections to the two-nucleon amplitudes 
or the reaction dynamics. The latter will introduce minor 
corrections to the diffraction-stripping term only, that 
contributes to the cross sections at the 30% level. Since 
the magnitude of these effects is not quantified here, 
we do not draw detailed conclusions regarding the ab- 
solute cross sections and Rs{2N) for ^^C, although the 
required corrections will not be large. We show the ra- 
tio (Texp/cfth in Table Hill to compare the results for the 
relative strengths of the different two-nucleon removal 
channels. All three channels will be affected similarly by 
the corrections mentioned above. 

We note that in the full calculations there is an en- 
hancement in the theoretical np removal cross sections 
for all three incident energies, already in excess of the 
16/6 ~ 2.7 ratio from the counting of pairs. The cal- 
culated np : pp cross section ratios, cr„p/crpp ~ 2.9, are 



TABLE IV: Like and unlike two-nucleon removal cross sec- 
tions (in mb) for a ^^C projectile incident on a carbon target 
at 2100 MeV per nucleon. The excitation energies, Ef, of 
each final state are shown in Fig. (2] The TNAs used were 
calculated using the wbp interaction. The sums show the ac- 
cumulated cross sections that lead to the ground state and the 
7-decaying bound excited states of the mass A= IQ residues. 



Residue 




T 


Ostr 


Ods 


OMf 


0-2N 


10 


0+ 


1 


1.59 


0.64 


0.06 


2.30 




2+ 


1 


1.96 


0.71 


0.06 


2.74 












sum 


5.04 












exp. 


4.11±0.22 


lOBe 


0+ 


1 


1.65 


0.68 


0.07 


2.40 




2+ 


1 


2.02 


0.74 


0.07 


2.83 




2+ 


1 


0.88 


0.32 


0.03 


1.23 




0+ 


1 


0.04 


0.01 


0.00 


0.06 












sum 


6.52 












exp. 


5.81±0.29 


lOg 


3+ 





5.11 


2.00 


0.20 


7.30 




1+ 





2.47 


1.01 


0.10 


3.58 




0+ 


1 


1.62 


0.66 


0.07 


2.35 




1+ 





1.81 


0.69 


0.07 


2.57 




2+ 





0.63 


0.24 


0.02 


0.89 




3 + a 





1.14 


0.43 


0.04 


1.62 




2+6 


1 


1.99 


0.72 


0.07 


2.33 




■y + a 





0.30 


0.10 


0.01 


0.41 




2+a 





0.75 


0.28 


0.03 


1.05 












sum 


19.02 












exp. 


35.10±3.40 



"states decay by o-emission with a 100% branching ratio, 
''state decays by a emission with a 16% o-branch. 



marginally smaller than those for the np : nn case, « 3.8, 
largely reflecting the different fractions of spectroscopic 
strength to bound states in the two like-nucleon cases - 
there being only two bound states in ^'^C. 

Since the magnitudes of the pp and nn removal cross 
sections are reasonably described, our expectation is that 
the cross sections to the states of the T = 1 isospin mul- 
tiplet in ^'^B are similarly well determined. However, an 
independent measurement of these cross sections would 
provide a useful verification of the direct nature of the 
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reaction. Thus, within the direct reaction model used, 
we must attribute the cross section deficit to the calcu- 
lated yields of the T = final states. We also note that, 
whilst the energy dependence of the like-pair removal 
cross sections is reasonably reproduced by the theoretical 
calculations, there is some discrepancy in the unlike-pair 
removal channel - with the result that the degree of un- 
derestimation varies with energy. The theoretical values 
for both the like- and unlike-pair removal decrease mod- 
erately with increasing projectile energy. This said, there 
remain considerable variations from the experimental np 
knockout cross sections. 



D. Momentum distributions 

In addition to the fragment production cross section 
measurements of Ref . [13| , residue momentum distribu- 
tions were also measured [36|. Except for the case of 
two-proton removal using a ^Be target, the experimental 
data sets have been published only as the widths of Gaus- 
sian fits to the experimental data. Here we denote these 
as Ag- Further, the published width parameters were 
the averages of data measured for several targets (Be, 
CH2, C, Al, Cu, Ag and Pb), though it was observed 
that, within the accuracy of the experiment, there was 
no dependence on the target mass above a 5% level. The 
references above do not discuss the thicknesses of the 
targets used, or of the momentum resolution of the in- 
cident beam, both of which could broaden the measured 
widths. We assume, since the data are for high-energy 
(stable) beams, that this experimental broadening of the 
width is negligible. 

Follo wing , and to enable comparisons with the results 
of Ref. [3al , we fit our calculated (projectile rest frame) 
inclusive residue momentum distributions with Gaussian 
functions. The two-nucleon stripping and diffr active- 
stripping mechanisms are expected to yield very simi- 
lar residue momentum distributions [17| . and here we 
present calculations for pure stripping events, scaled to 
match the mechanism-inclusive cross section for each fi- 
nal state. In general these calculated final-state-inclusive 
residue momentum distributions are very close to Gaus- 
sian shapes. This is not the case however for the individ- 
ual (exclusive) residue final states, particularly those of 
higher spin which have a broader and fiatter distribution 
near the central value, Kc — MeV/c. The results are 
collected and compared in Table |V] and are, encourag- 
ingly, in good agreement. 

The calculated residue momentum distributions for 
two-proton removal are compared to the data in Fig. 
Ip the data having been scanned from Figure 1 of Ref. 
[36| . The data provided are for a beryllium target while 
the calculations assume a carbon target. Any differences 
in the momentum distributions are expected to be very 
small for these two light targets. Both the exclusive con- 
tributions from the four ^"^Be final states involved and 
their sum are shown. The figure shows clearly that the 



TABLE V: Projectile rest frame residue momentum distri- 
bution widths Ag (for reactions at 2100 MeV per nucleon) 
obtained by fitting Gaussian profiles to the experimental 3y] 
and theoretical (this work) inclusive distributions to bound 
final states. The TNAs used were calculated using the wbp 
interaction. The data of Ref. [s^ are the average widths from 
data for several targets (see text). The calculated values are 
for a carbon target only. 



Residue 


A ea-p 


Ath 


"B 


106±4 


99 


"C 


103±4 


100 


lORe 


129±4 


127 


lOg 


134±3 


132 




121±6 


120 




(MeV/c) 



FIG. 4: (Color online) Comparison of measured and calcu- 
lated ^''Be projectile frame residue momentum distributions. 
The red curves show the ^°Be final-state-exclusive results for 
the ground state (solid line), 2f (dashed line), 23" (dot-dashed 
line) and Oj (dotted line). The inclusive cross section, their 
sum, is shown by the black solid line. The data is for a ^Be 
target whereas the calculations assume a ^^C target. The cal- 
culations have been offset by —30 MeV/c and the have been 
scaled to match the experimental two-proton removal cross 
section on the ®Be target, of 5.97 mb. The TNAs used were 
calculated using the wbp interaction. 



inclusive Gaussian-like distribution is formed of exclu- 
sive cross sections of very different widths for the differ- 
ent final states. Thus, residue momentum distributions 
for these individual final states would be of considerable 
value in validating, or revealing deficiencies in the de- 
tails of the shell-model wave functions and TNAs used 
here, including our description of the reaction as a direct 
single-step process. 

Details of these calculated full width at half maximum 
(FWHM) widths in the case of np-removal at 2100 MeV 
per nucleon, for the six 7-decaying final states of ^°B, are 
shown in Fig. [S] Results from both the wbp (solid lines 
and filled points) and pjt (dashed lines and open points) 
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FIG. 5: Calculated exclusive cross sections (a) and full width 
at half maximum (FWHM) widths of the momentum dis- 
tributions (b) for the six 7-decaying final states of the '^''B 
residues, for reactions at 2100 MeV per nucleon. The TNAs 
used were obtained using the wbp interaction (solid lines 
and solid points) and the pjt interaction (dashed lines and 
open points). The significant variations predicted: (i) for the 
widths of the momentum distributions to the different final 
states, and (ii) in the sensitivity of the ^°B, T = final state 
yields to the effective interaction used, are evident. 

shell-model interaction TNA are shown. The widths from 
the two interactions are remarkably similar for all transi- 
tions, indicative that the XiS-contributions to each final 
state are very similar [l^ . The momentum distributions 
for the first and second T = 0, 1+ states of ^"B were also 
presented in Figure 3 of Ref. [l^, calculated from the 
WBP interaction. The corresponding exclusive cross sec- 
tions (upper panel of the figure) also show very similar 
trends for the two shell-model interactions. They differ 
in detail, however, particularly with regard the ratio of 
the cross sections to the ground (3+,r = 0) and first 
excited (1^,T = 0) states. We note also that the two 
interactions predict very similar cross sections in the two 
T = 1 states, the differences being more evident in the 
T — states yields. Once again, this suggests interesting 
information would be gained from exclusive final states 
yields. 

The predicted final state-dependent variations of the 
FWHM widths from the model used also offer consid- 
erable scope for further testing of the reaction inputs, 
the shell-model wave functions used, and of the nucleon- 
nucleon correlations content of these wave functions. 



V. SUMMARY 

We have considered in detail theoretical expectations 
for the cross sections of two-nucleon removal reactions, 
from ^^C incident on a carbon target, at beam energies 
of 250, 1050 and 2100 MeV per nucleon. The calculated 
inclusive cross sections for two like-nucleon (T = 1) re- 
moval are broadly consistent with available experimental 
data. For np-pair removal the analogous T = 1, 
final states are also shown to be insensitive to the shell- 
model interactions used. The calculated inclusive np-paii 
removal cross sections at all energies underestimate the 
data by approximately a factor of two. Theoretical cal- 
culations of the widths of the final-state- inclusive residue 
momentum distributions on the other hand are consis- 
tent with the available experimental data, including for 
the np-removal channel. 

Further measurements, of final-state-exclusive cross 
sections and residue momentum distributions, would al- 
low a much more detailed scrutiny and validation of the 
direct nature of the reaction, including the identification 
of any indirect reaction components. The calculated np- 
removal cross sections to the T — 0, final states were 
also shown to have sensitivity to the shell-model effective 
interaction used; for example, the ratio of the calculated 
cross sections to the ^"B ground (3+,T — 0) and first 
(1+,T = 0) excited states. 

The overall conclusion from the present analysis is that 
the existing data suggest that the T = 0, np-spatial corre- 
lations present in the wave functions used are insufficient. 
We have shown that exclusive measurements would offer 
a means to interrogate these shell-model inputs, in partic- 
ular the np-channel, T = wave functions, and the direct 
reaction mechanism predictions in considerable detail. 

A similar two-nucleon knockout study can also be per- 
formed using an ^^O projectile, for which np-pair knock- 
out using electromagnetic probes has been measured 
[13, and that supports the SRC observations for ^^C. 
In this case the ^'^O residue is also of interest, given that 
only the 0+ g round state is bound, and where the results 
of Refs. |l3l 136} indicate a correspondingly small cross 
section and a narrow momentum distribution, consistent 
with theoretical expectations. Additionally, possible 2hLo 
components in the ^^O ground state wave function add 
interest to the magnitudes of the absolute cross sections 
in this case. 
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